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I. Introduction : 

One considers a 2-dimensional model for an electron of electric charge e and of mass 
m subject to a random magnetic field. Here, random magnetic field means a Poissonian^ 
random distribution {r^, z = 1,2,...., N} of fixed infinitely thin vortices carrying a flux (p, 
modeling magnetic impurities, and characterized by the dimensionless Aharonov-Bohm 
{A — B) coupling a = e0/27r = 4)/(l)o- 

The question we ask is about the effect of disorder on the energy level density of 
an electron -test particle- averaged over the random position of the vortices Q. In the 
thermodynamic limit N oo, y — > oo for a distribution of vortices of density p = N/V, 
one might naively argue that the average magnetic field, < B >= apcpo, is expected to 
become meaningful in the limit p — > oo, a ^ 0, with pa kept finite. On the other hand, 
if p is finite, and a non vanishing, corrections due to disorder should exhibit non trivial 
magnetic impurity signatures, like broadening of Landau levels and localization. 

In a first paper ||2| on this problem, two approaches were used: 

- A path integral Brownian motion analysis where the problem at hand is mapped, 
after averaging on disorder, to a study of winding properties of Brownian curves. 

- A quantum mechanical formulation, where the contact hard-core boundary conditions 
at the location of the impurities where properly taken into account by an appropriate wave 
function redefinition, allowing for an analytical averaging on the random singular A — B 
interactions. 

The main results were : 

- in the presence of a constant external magnetic field, if the total (external -|- average) 

magnetic field is strong enough so that one can neglect the coupling between the low- 
One could consider as well other types of probability distribution - Gaussian for example. For the 
sake of simplicity, we concentrate only on the Poisson case. 
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est Landau level and the excited Landau levels by the random component of the vortex 
distribution, the system was best described when projected in the LLL. Since one has in 
view a sufficiently dilute gas of electrons compared to the available quantum states in the 
LLL -the fractional Quantum Hall regime-, such a restriction is licit. In this situation, the 
quantum problem was explicitly mapped on a problem of random 6 impurities, where the 
average density of states happened to be known 

- in the average magnetic field limit a ^ 0, p ^ oo, a global zero point energy shift 
^^^^^ =< uJc > materialized in the Landau spectrum of the average magnetic field < B >. 
The origin of this shift was traced back to the hard-core boundary conditions at the 
location of the impurities |^. 

- in the case a = 1/2 and p finite, on the other hand, Brownian numerical simulations 
showed that the sole effect of the impurities was a depletion of states at the bottom of the 
spectrum, characterized by a Lifschitz tail in the average density of states, caused by the 
isolated impurities. 

An additional interesting result was that the average density of states < p{E, a, p) > 
happened to be a function of E/p and a only (this is also true in presence of an external 
magnetic field). The < p{E/p,a) > scaling implies that the impurity density p is not a 
relevant parameter, and can be arbitrarily set to a given value, since changing its value 
amounts simply to a rescaling of the energy unit. A consequence is a more precise definition 
of the limit where the average magnetic field becomes meaningful. It should happen, 
for any given finite p, when a becomes sufficiently small, independently of 
the value of p. Thus one expects a critical value Oc where a transition occurs 
from an almost free density of states with a Lifschitz tail at the bottom of 
the spectrum (oc < a < 1/2), to a Landau like density of states, with Landau 
oscillations, i.e. Landau levels separated by a Landau gap (0 < a < Oc). 



3 



A semiclassical understanding of this transition consists in taking an electron with a 
Fermi velocity vp, subject to the average magnetic field < B >= pa(f)„. A typical cyclotron 
orbit radius is i? ~ vp/ipa)- For a given p, one has ~ 1/a, thus the smaller is a, the 
larger is the number of magnetic impurities enclosed by the cyclotron orbit, and therefore 
the more accurate is description in terms of the average magnetic field. To summarize, 
for a given finite p, the smaller is a, the better is < 5 >, eventhough it is smaller and 
smaller. However, there is no other magnetic scale in the problem to compare it with. 

An important consequence of the Lifschitz tail - Landau levels transition should be a 
non conducting-conducting transition, since one expects localization due to disorder when 
a > ac, and extended states in the Landau regime when a < ac- 

The aim of the present paper is to show that the critical value at which the transition 
occurs is ac — 0.3 — 0.4. 

We will present two types of evidence for this result : 

- Numerical evidences 

i) for the specific heat, where a Brownian motion numerical study indicates a transition 
from an almost free type to a Landau type specific heat when q""™ ~ 0.28. 

ii) for the density of state, where, under some reasonnable simplifying assumptions, 
the transition is shown to occur at ac — 0.35. 

- Analytical evidences : 

i) one can explicitly show that the specific heat transition is possible only when Landau 
oscillations have already begun, implying that ac> ag, but close to ao- 

ii) quantum mechanical evidence using the impurity cluster expansion of the average 
partition function for an arbitrary large number of impurities. At a given order p^ , for a 
test particle subject to a given number N of impurities described by the Hamiltonian iJ^r, 
the average partition function < Tr exp(— /3i7jv) — Tr exp{—(3HQ) > has been computed 
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as a power series in a. One recovers, at order the leading order in a, the partition 

function of the average magnetic field. Higher order corrections can be computed : for 
N = 2, the diagrammatic expansion has been done up to order a^. We will show that 
this is all what is needed to obtain ao — 0.29, which is indeed very close to the numerical 
result a™"" = 0.28. 

II. Definition of the model : 

The system is peridodic in a with period 1, so one can always take a G [0,1]. In 
the absence of an external magnetic field, there is no privileged orientation to the plane, 
therefore the system is symetric wrt a = 1/2, and one can restrict a S [0, 1/2]. It follows 
that any physical quantity of interest should depend on a(l — a) only. 

II. a Path integral Brownian approach : 

One starts from a square lattice of N squares of size a^, in which point magnetic 
impurities are randomly dropped. Let Ni be the number of vortices dropped on square i. 
A random configuration {Ni} will be realized with the probability 



with N/M = Ylh=i ^i/-^ = In the thermodynamic limit, this is nothing but a Pois- 
sonian distribution. In order to compute the average level density < p{E) >, one focuses, 
in the thermodynamic limit N ^ oo, on the one-electron average partition function per 
unit volume (for an electron of unit mass and charge) 



where {C} is the set of L steps closed random walks, and rii is the number of times the 
square i has been wound around by a given random walk in {C}, i.e. its winding number. 
Zo = 2^ is the free partition function per unit volume, with t (the inverse temperature 13) 




(1) 



Z = Zo< e*E.=i2^">^- 



>{C,Ni} 



(2) 
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the length of the curve (2t = La^, e = m = 1). (||) is obviously invariant when a is shifted 
by an integer, and when a — > —a (because rii comes always with —rii), so one can always 
restrict to0<a<l/2. Averaging Z with (||) one gets 

Z = Z,<e^i:.^"(^"--i)>|C} (3) 

where Sn stands for the arithmetic area of the n-winding sector of a given random 
walk in {C}. Eq. (^) is still true in the limit a ^ 0, L — > oo, with t fixed, i.e. for 
Brownian curves in the plane, yielding a path integral formulation for the problem at 
hand. Extracting the variable t, it rewrites as 

Z = ZoJ e-P^'^^+'^^PiS, A)dSdA = Zo< e-^*(^+^^) >{c} (4) 

where 5 and A are defined as 

2 

S = - ^ Sn sin (vran); < S >= 7ra(l — q) (5) 
* n 

A = -J2Sn sin(27ran); < yl >= (6) 

* n 

P{S, A) is the joint probability distribution of the random variables A and S. One has 
used the general property of Brownian curves that 5„ scales like t (remind that < Sn >= 
t/{2'Kn?) 1^]), so that the variables S and A are actually t independent. This implies that 
the average partition function Z in has the form F{pt)/t, and thus its inverse Laplace 
transform, the average density of states, is necessarily a function of E/p and a. 
II. b Microscopic Quantum Hamiltonian : 

The Hamiltonian of an electron of mass m and charge e subject to the potential vector 
of N vortices at position fi (k is the unit vector perpendicular to the plane) reads 0] 



„ 1 / ^ A; X {r-rij ^ 
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Because of periodicity in a, one can always restrict a to a G [—1/2, +1/2]. We consider 
in the thermodynamic limit the Poisson probability distribution dP{fi) = dfi/V. 

The system described by the Hamiltonian (0) is not yet entirely defined. Boundary 
conditions on the wave functions have to be specified when the electron comes close to an 
impurity. This can be achieved in a non ambiguous way by adding to the pure Aharonov- 
Bohm Hamiltonian the contact terms 

-*^i(-i:"%5p)^±i:^^v-n) (8) 

These contact terms amount to couple the spin-up (+) or spin-down (— ) degree of freedom 
Q of the electron endowed with a magnetic moment /i = — 2m 1*^1 (thus an electron 
with gyromagnetic factor g = 2), to the infinite magnetic field inside the flux-tubes. 
Their origin can also be understood in a perturbative framework : the original Aharonov- 
Bohm spectrum with vanishing wavefunction at the location of the magnetic impurities 
(a particular self adjoint extension describing impenetrable vortices : hard-core) can be 
perturbatively obtained if and only if the contact terms with the (+) sign are taken into 
account, whereas the contact terms with the opposite (— ) sign correspond to a different 
self-adjoint extension, where the the wavefunction is singular at the location of the vortices 
(the particle is attracted inside the vortices : attractive-core). Note that in the Brownian 
path integral formulation, only the hard-core case can be described, due to the fact that, 
by definition, a given Brownian path has no chance to pass through a given impurity 
location. 

The contact terms happen to be crucial for the averaging on the disorder : consider 
the nonunitary wavefunction redefinition Q 

N 

i=l 
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to obtain the Hamiltonian acting on V'jv(^) where the impurity potential now reads 

^^=--9.a.-^f:^+^^f:^ (10) 

ra m ^ z — Zi m ^ z — Zj 

1=1 ' 1=1 * 

which rewrites if a > as 

2 „ „ 2a 



H^ = d,d,-—Y,^^ (11) 

m m ^ z — Zj 

1=1 ' 

^^=--5.a.+^f:^ (12) 



and if a < 



m m " z — Zj 

1=1 



2 „ ^ 2aS^ d. 



H+ = d,d, + — Y,^ (13) 



m m ^ z — Zi 



2 _ 2a A a 



H-^ = d,d,-—Y^^^ (14) 

m m z — Zi 

One sees in both Hj^ cases that there is again no distinction to be made between a £ 
[0, 1/2] and a E [—1/2,0]. Indeed, the impurity potential are hermitian conjugate the one 
from the other, but perturbative computations in a are always such that real results are 
obtained. So one can always restrict to a € [0, 1/2]. Both cases should be considered 
in principle. However, it is easy to see that they are equivalent, except for a global shift 
+ < > in the hard-core case, and — < Uc > in the attractive-core case. 
Indeed, consider, instead of (^, the wavefunction redefinition where the average magnetic 
field Landau pre-exponential factor has been extracted Q 

N 



V'^(f) = e^2-<-^>- Yl \r- r-;|±°^±(f) (15) 

i=l 

to get (in the thermodynamic limit ^ oo) 

H^ = ±<u;c> +H±<B> + V^{a)- < V^{a) > (16) 



where H^b> is the Landau Hamiltonian for the average < B > field and the impurity 
potential reads 

2a ^ d ^ z 

V^(a)- < V^(a) >= ^_ +2 <u;c> zdz+Y^ <u)c> a- — = — m < Uc >^ zz 

mf^^z-Zi ^ z-Zi 

(17) 

{V~{a)— < V~{a) > is obtained by taking the hermitian conjugate of V~^{a)— < 
y+(a) > and a — —a). Clearly, in the limit with no disorder a — 0, one obtains 
the Landau Hamiltonian for the average < B > field with a ± < Uc > shift. 

The Hamiltonians and are equivalent, and can be indifferently used for com- 
puting the partition function or the density of states. However, interactions with two 
magnetic impurities have disappeared from H^, greatly simplifying the average on the 
disorder, which can be easily done using the identities 



/ 



/ 



/ 



dzidzi- — —dz = TTzdz (18) 

Z- Zi 

11 z z' 

dzidzi- — -dz- -dz' = Tr{- z + - — z;)^^^^' (19) 



Z Zj^ Z Zi z z z z 

dzidzi _ ^ _ dz dz' _„ \ _ dz" 

Z Z-i Z Z-i Z Zi^ 

J 



''\z' -z)(z^^ -z)^ iz^^ --z<){z--z<) + (z'-z")(z--z-")^^^^^'^^" ^^^^ 

etc.... 

Simple dimensional arguments at the level of the Hamiltonian are sufficient to under- 
stand the scaling property of the average density of states. Rescaling the length unit by 
A amounts to rescale the Hamiltonian, thus the energy, by 1/A^. On the other hand, in 
d = 2, the same length unit rescaling implies for the density p the same 1/A^ rescaling. It 
is not surprising to find, after averaging, the Ej p scaling of the density of states. 

In the sequel one will concentrate only on the hard-core case, bearing in mind that the 
attractive-core case can be straightforwardly deduced from the hard-core case. 



II. c Known results: the cases a and a = 1/2 

To simplify notations, we will use from now on the Brownian notations e = m = 1, 
and (3 = t. 

From IQ] we know that : 

i) In the limit a ^ 0, one expects from @ that Z Z^b> = Zo < e*<'^>I^n "'S'" >^^^^ 
the partition function of one electron in an uniform magnetic field < B >= apcpo- However, 
possible corrections coming from the exponent exp(i27rna) — 1 might alter this result. 
Due to the non-differentiability of Brownian paths, X^n^^'^'n is not defined for a typical 
Brownian curve where < Sn >= t/(27rn^) |p. Certain recent results in the mathematical 
litterature Q show that, for n sufficiently large, n'^Sn — >< n'^Sn >= t/(27r). It follows that, 
when a^O, 5'~<5'>~ ira, and A ~ 2TraA, where A is the algebraic area enclosed by 
the Brownian curve {A = '^nSn/t). One deduces that Z — >q,_^o ■^<B>e~*^^^/^, implying 
that the system of random vortices is equivalent to an uniform magnetic field < B >, but 
with an additional positive shift in the Landau spectrum < B > /2 =< lOc >■ Note that 
Z<B> is built by the random variable A, and the shiff by the random variable S. 

ii) When a = 1/2 on the other hand, one can explicitly test the effect of the random 
distribution of vortices. (^) now reads Z = Zq < e~^^^ ^{C} ^^^d the average density of 
states, obtained by inverse Laplace transform of Z , is 

r- 

< p{E) >= po{E) / " P{S)dS (21) 

P{S), the probability distribution for the random variable S, was estimated numerically 
by simulations on a lattice, where a number of steps ranging from 2000 to 32000 was used. 
In Figure 1, < p{E) > displays a Lifschitz tail at the bottom of the spectrum, around 
E p < S >= vrp/4, where a behavior < p{E) >~ exp(— p/£') is expected. 

In both these extreme cases a = 0, a = 1/2, an energy level depletion at the bottom of 
the spectrum is observed (global shift in the Landau spectrum in the former case, Lifschitz 
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tail in the latter) . This pattern is quite reminiscent of one impurity A — B density of states 
depletion pi{E) - po{E) = ^^^^6{E) @, which occurs precisely at E = Q. 

The question we now ask is what happens when a continously decreases from a = 1/2 
to a = 0? In particular, is it possible to understand the transition from a Lifschitz tail 
pattern to a Landau pattern in terms of several impurities average density of states? 

III. Numerical evidences 

III. a Specific heat 

Interesting enough is the study of the specific heat averaged over disorder 

c = ^*'^<l^^'>W} (22) 
where k is the Boltzman constant and 

Z' = Zo<e'^ >{C} (23) 

is the partition function for an electron in a given distribution of vortices {Ni}. We first 
note that 

< InZ' >{^j= ln[< Z' >{ivj](= InZ) (24) 

In principle this property holds only for short ranged impurity potentials. This is of 
course not the case in the present problem. However, before averaging over disorder, the 
partition function involves only closed Brownian curves, and thus is entirely determined by 
the impurities distributed inside each Brownian curve (see @). Using the basic property 
that two Brownian curves on the plane have no chance to intersect each other, it is easy 
to get 

< Z'2 >|^,}=< Z' >\^^^ (25) 
i.e. Z' is a self- averaging quantity. A straightforward computation leads to the t (inverse 
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temperature) expansion (note that P{S, A) = P{S, —A) 



c ~ Co + kt\< >{c} -<S >lc} -<A^ >{C}) + ■■■ (26) 

Both quantities < >{C} — < S >^c} ^^'^ ^{<^} '^^^^ ^ natural interpretation 

in terms of a Landau spectrum : < A^ >{C} determines the average magnetic field Landau 
levels, and < >{C} ~ < S >f^c} measures the deviation from the Landau pattern due to 
disorder. It is not a surprise that their difference (< S'^ >{C} — < "S" >{c})~ -^{C} 
plays a role in finding the critical point for the specific heat. They have been studied 
numerically, and are displayed in Fig. 2 for 2000 random walks of length L = 100000. 
When a continously decreases from a = 0.5 to 0, c — Cq is first positive, then vanishes 
for a ~ 0.28, then becomes negative. The vanishing of the t'^ term in the specific heat 
corresponds to a perfect gas behaviour. We will argue below, by analytical means, that 
the perfect gas behaviour takes place at — 0.29. 

Ill.b The density of states : 

In order to estimate < p{E) > numerically, a simple heuristic ansatz has to be made. 
Correlators between the variables S,A and the algebraic area A = J2n''^'^n/i, such as 

_ < SA> - < S >< A> 

^'^ " ^/i<S^> -<S >^)(< A^>-<A>^ ^ ^ 

or Cs,A a^re trivially vanishing because of obvious symmetry properties, and carry no in- 
formation. It is more appropriate to use the variables \A\, and |^|. Numerical simulations 
indicate that (^^1,1^1 — 1 for a < 0.05, and (^1^1,1^1 — 0.93 — 0.95 for a ~ 0.25. One 

deduces that the variables \A\ and |^| are highly correlated, in particular when a is very 
small, i.e. in the average magnetic field limit. We thus assume the linear relation 

1^1 = (28) 
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where B is an effective magnetic field determined hy B/p = y ^^^^ = Vl2 < >, which 
gives, in the hmit a — > that B = 2Trap =< B >, as it should^. On the other hand, in 
the limit a ^ 1/2, B = 0, the magnetic field is absent, as again we already know from 
(^. In conclusion, the ansatz ( |28D seems quite reasonable. Now, introduce the positive 
random variable 



S' = S- Cs,iAnl 1,1^. 1-^1 (30) 



< g2 > _ < g >2 

<A^>-<\A\> 

so that S' and |^| are uncorrelated. The average partition function rewrites as 



Z^Zo< e-P*^' >< ef^^^'-^^ cos{tBA) > (31) 
The inverse Laplace transform of Z/ < e^^*^' > can be obtained analytically 



1 - 

1 + e-2" + 2e-" cos 7 



P(E') = Po{E') ^ , , , (32) 



where 



-,=2.mmA (33) 



and 

7/x 



(5/p) 



(34) 



where p = C'gj^i of^^>'-<j^f>^ • "^^^ average density of states rewrites 

cE/p 



2 



< p{E) >= / P{S')p'{E - pS')dS' (35) 

JO 



''To obtain < A^ >= 1/12, consider simply the probability distribution for the random variable .4, 
which is nothing but the Fourier transform of the partition function of a charged particle in a constant 
magnetic field |^ 



1 ,„ ,R4 B 



P(A) = — / dSe'-^-^ — (29) 

^ ' 27r / 2sinhB/2 ^ ' 

J —00 / 
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where P{S') is the probabihty distribution for the variable S' . The numerical results are 
displayed in Fig. 3 for different values of a. It is clear that the average density of states 
is no more monotonic when a is below the critical value ac — 0.35. Also, when a — ^ 0, 
well separated Landau peaks appear with a well defined Landau gap. In this low energy 
region (which corresponds to long Brownian curves), the system clearly mimicks a constant 
average magnetic field. 

IV. Analytical evidences 

IV. a Oc > Oo : 

Let us show that the transition for the density of states occurs necessarily at a critical 
ac > Q.O- Let us assume that when E — > oo, < p{E) >— > Po{E) = i.e. at very high 
energy the system does not see the magnetic impurities, and when E ^ < p{E) >— 0, 
i.e. at very small energy the depletion of states due to the hard-core impurities is effective. 
Both these assumptions are quite reasonnable, and are actually verified in all the numerical 
simulations and analytical studies. After integration by part in Z, one gets 

, kt' /o°° /o" d^d£;-e-^(^+^')^<^^ -^<^(g)/^> (E - E'f 

Uo ^ dE i 

At small i, this expression becomes 

c . c. + ktv iEiE''t^£mL±^^lJim!l2i^E _ ef (37) 

which in turn implies that for c — Cq to change its sign one necessarily should already have 
Landau oscillations. 

IV. b The concentration expansion : 

In the thermodynamic limit, the number of impurities is only fixed on average. We are 
in a situation similar to the one encountered in statistical mechanics for an undetermined 
number of particles, where there exists a cluster expansion for the thermodynamical po- 
tential, which involves at a given order of the fugacity , the partition functions for N 
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particles, — 1 particles, etc... Here, there is one particle subject to the effect of an 
undetermined number of impurities. Again, there exists an expansion for the averaged 
one particle partition function [|lO|, analog to the cluster expansion. The parameter of 
the expansion is the density (concentration) of impurities p, and at a given order , the 
partition function per unit volume is expressed in terms of the partition functions for 
impurities, — 1 impurities, etc... It reads 

Z = Z,{l + pV{^-l) + y^V^[^-2^ + l] + y^V\^-^^ + ^^-l] 

Zjq l\ Zjq Zq 6\ Zq Zq Zq 

where =< Tr exp{—tHj^) > /V. 

Let us compare (^) with the pt expansion of (4), namely 

Z = Zoil-pt<S> <iS + lAf > -^f^t^ <iS + iA)^ >+•••) (39) 

So, considering p as the expansion parameter, one sees that in order for both expansions 
(^) and ( |39|) to match one necessarily has 

|-1 = -A<5>. |-l-2(|-l) = |l<(5 + .Af>.... (40) 

One already knows that 

i) when a ^ 0, (^) is the partition function per unit volume of a charged particle in 
the average magnetic field, with a Landau spectrum shited by < Uc > 

1 <B> ^-t<B>/2 1 ^C(-n) 



-e 



-t<B>/2 ^ _:_(-^_sr^ ^\ ^B>tr+^) (41) 

' n=0 

(again this follows directly from general properties of Brownian curves -note also that the 
attractive-core case would be obtained simply by<5>— > — <-B>in (^)). 
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ii) from the t expansion of the specific heat c = Cot^ '^ ^^i'^ that the Oo value is attained 
when 

< {S + iAf >=< S (42) 

It follows that by a diagrammatic expansion in a of quantities such as ^ — 1, ^ — 
1 — — 1), • • •, one should recover 

i) the leading a behavior ( p| ) 

ii) the critical ao value by finding the zeroes of(^ — 1)^ = ^ — 1 — 2(^ — 1) 

The diagrammatic expansion of Z]\j is performed by usual perturbative methods for the 
Hamiltonian (|ll|). One uses a non perturbed basis of free thermal propagators Gt{r,r') = 
2^ exp( — (r — 7^)^/21). In the computation of the average impurity partition function, 
one encounters, at a given order in a, volume divergences from the dzidzidzjdzjd... space 
integrals over the impurity locations fijfj,- ■ if computed directly in the thermodynamic 
limit. Thus, in principle, the need of a certain regularization prescription, as for example 
a harmonic regularization. However, in our case, one can circumvent this difficulty by 
computing, rather than average partition functions Tr exp{—tH]\f), the average of the 
thermal propagators G]\f{ro,ro) =< fo\exp{—tH]\f)\fo > from and to a given point Tq. 
The actual average partition functions are by definition space integral over Tq of these 
propagators. Yet, averaging over disorder has to be made. If it is done before the final 
fo space integration, it follows that < GAr(ro,ro) > does not depend on fo anymore. 
Therefore, the final space integration becomes trivial, since it amounts to multiply by V 
the coinciding point propagator. In this computational scheme, which relies on the crucial 
fact that averaging on the disorder has to be made, there is no need for a particular 
regularization procedure, since the infinite volume factors out trivially in the last space 
integral. 

At a given order a", a diagram has n vertex, and m < n impurities: it contributes at 
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order ot^ in the concentration expansion. The diagrams which contribute to the < B > 
average partition function are necessarily such that m = n, and corrections to the average 
field limit are built by diagrams with m < n. We have analytically computed the leading 
< B > diagrams up to order a^/j^, and at order p'^ -the two impurity case-, the corrections 
to the leading a'^p'^ digram up to order a^. The diagrammatic expansion results are: 



1 a{a — 1) 
V%, 2 
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^-l-2(^ 
Zq Zq 

where ^ - ^C{3) = -0.48423323, 



1 



,1 



1 



7 



24 16 



C(3))a^ + 



Zs 
Zo 



Zo z^ 

Zq Zq 



{VZo) 



r(Oa^ + ---; 



(43) 
(44) 

(45) 



Za , Z-i Zi , Z^ 1 

— - 4— + 6— - 4— + 1 = -I 

Zn Zn Zn Zn iyZ^Y 



1 



0^ + 



- . . (46) 

1/0 z/o z/o \\' ^o)^ 30 

We see that at leading order a^p^, (43, 45, ^6|) indeed reproduce the average magnetic 



field expansion (41). 

We are now in position to determine the value at which the specific heat transition 
occurs, just by considering (43,4^ on the one hand, (4^,42|) on the other hand. Before 



doing so, the diagrammatic expansion (44) should be first completed to obtain a a(l — a) 
polynomial, since this should be so. One gets 

1 fl 



Zo 



Zo ' -^Zo (y^j2vi["(l-«)]'+^["(l-«)]'+^(l-^^(3))[a(l-a)]^+0([a(l-a)]5) 



1 , 1 
2' 



The resulting polynomial equation in a writes 



a(l — a) 



1 



1, 



-[a(l - a)Y + -[a{l - a)Y + -(1 - -C(3))[a(l - a)]' 



(47) 



(48) 
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In the interval a G [0, 1/2], this equation has the desired root Uq = 0.29. 

One can go a bit further, and compare the diagrammatic expansion with the numerical 
estimations for (< 5^ >{C} — < S >{c})~ < ^ ^{<^})' -^^S' ^' "-"^^^ three 
curves - a) and b) are numerical simulations, c) is the diagrammatic expansion- have the 
same intercept with the horizontal axis (thus the same Uq ~ a""*" value), but also they 
exhibit the same qualitative behavior. Differences are due to the fact that actual winding 
properties for true Brownian walks are difficult to reach from finite length random walks, 
and to the fact that the perturbative expansion in a(l — a) is by definition incomplete. 
Still, it is remarkable that both approaches yield the same value for ao, altogether with 
the same a qualitative behavior. A convergence between both approaches is expected. 

An interesting consequence of the discussion above concerns the average density of 
states for the one impurity problem, two impurities problem, etc. In the simple one 
impurity case, averaging over disorder is trivial, since the Aharonov-Bohm low energy 
depletion of states does not depend on the position of the vortex. One simply recovers the 
standard result Q 

< p^{E) - po{E) >= p,{E) - po{E) = ^^^^5{E) (49) 

In the two-impurities case, things become highly non trivial, but averaging has allowed 
for the diagrammatic expansion result (|4^). One deduces that 

< p^{E) - po{E) > - 2 <p,{E)-po{E)>='^{^^[a{l-a)]''+^-[a{l-a)f 

+ I(l-ic(3))Kl-a)r + ...)y(S) (50) 

<p^{E)-po{E) >-2,<p2{E)-po{E) >+3<pi{E)-po{E) >= iyf{0a^ + ---)5 (E) 

(51) 
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< Pi{E) - po{E) > -A < p-i{E) - po{E) > + 6 < p2{E) - po{E) > -A < pi{E) - po{E) > 

etc... At leading order p^a^, summing up all leading contributions from the one, two, • • ■ 
impurity cases, leads to 

V < B > °° 

= ^ <5(ii; - n < B >) - poiE) = pl{E) - po{E)i53) 

n=l 

which reproduces, as expected, the shifted Landau density of states Pl{E) — Po{E). 
IV. c Starting from the Landau basis: 

Perturbative computations can also be done by starting directly from the Hamiltonian 
(p!l|). The expansion of the average partition function will be modified in the following 
way: 

- propagators to be used are Landau propagators in the average magnetic field 

^ ,,, \ ^ — ^— \\z — z'Pcosht <iOc> +sinht < cJc > {zz — zz' 

' 27r sinh t < uJc> 

instead of free propagators. 

- the vertex to be used is 

-2a _ {dz -\<ujc>z) + 2a^{zd^ - <ujc> zz) 
z — z\ 2 N 2 

instead of 



-2a- 



Z - Zi 
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One has already noted that any diagram with isolated impurities only, i.e. of the a^p^ 
type, contributes to the average magnetic field partition function. It follows that, by defi- 
nition, it vanishes in the present formulation since the average magnetic field contributions 
are ab initio incorporated in the Landau propagator. It is possible to show that for any 
subleading diagrams a"/?'", m < n, the vertex reduces to 

-2a _^ _ {dz - I- <ujc> z) 
z — Zi 2 

Any such subleading diagram automatically incorporates the corresponding diagram 
in the free propagator approach, plus all diagrams deduced from it by adding an arbitrary 
number of isolated impurity lines, i.e. diagrams of the type a"^ p"^ (1 + ap+ a'^ p'^ + •••), m < 
n. 

V. Conclusion 

A more precise value of Oc and a better understanding of what is exactly happening 
at the transition is still missing. How Landau levels actually appear when a continuously 
decreases from Oc to 0? What is the exact nature of the transition? 

Another important issue concerns the conductivity and localisation properties of the 
test particle in the random system of magnetic impurities. In particular, one would like 
to have information on the Hall conductance in the impurity system. This question could 
be partially adressed by considering the average persistent currents due to the vortex 



distribution 1 11 1 . One has been able to get perturbative expansion in a for the one, two, • • ■ 
impurity average densities of states. In principle, one can deduce [11| from these density 
of states the average persistent currents for one vortex, two vortices, • • •. In the limit 
Q ^ 0, at leading order a"p", diamagnetic persistent currents for a constant magnetic 
field should be recovered by summing up all contributions coming from the impurities 
persistent currents. 
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Figure captions : 

Figure 1: The average level density of states < p{E) > as a function of the variable 
E/p for a = 1/2 exhibits a Lifshitz tail at the bottom of the spectrum. 

Figure 2: The variances of the random variables S and A obtained by simulations (2000 
closed random walks of 100000 steps) are plotted as a function of a. The intersection of 
the two curves at a"""* = 0.28 determines a change of the specific heat behavior at high 
temperature. 

Figure 3: The average level density of states, determined by numerical simulations, as 
a function of the variable E/p for different a values. When a < — 0.35, the density of 
states is no more monotonic and oscillations appear. 

Figure 4: The difference between the variances of S and ^ as a function of a, deter- 
mined by numerical simulations for a) 3000 closed random walks of 400 steps, b) 2000 
closed random walks pf 100000 steps, and also by the diagrammatic expansion to fourth 
order in a. The three curves intercept the horizontal axis at a ~ 0.28 — 0.30. The origin 
of the quantitative differences are explained in the text. 
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